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Generalized Bose-Einstein condensation (GBEC) involves condensates appearing simultaneously 
in multiple states. We review examples of the three types in an ideal Bose gas with different geome- 
tries. In Type I there is a discrete number of quantum states each having macroscopic occupation; 
Type II has condensation into a continuous band of states, with each state having macroscopic 
occupation; in Type III each state is microscopically occupied while the entire condensate band 
is macroscopically occupied. We begin by discussing Type I or "normal" BEC into a single state 
for an isotropic harmonic oscillator potential. Other geometries and external potentials are then 
considered: the "channel" potential (harmonic in one dimension and hard- wall in the other), which 
displays Type II, the "cigar trap" (anisotropic harmonic potential), and the "Casimir prism" (an 
Q . elongated box), the latter two having Type III condensations. General box geometries are consid- 

■ ered in an appendix. We particularly focus on the cigar trap, which Van Druten and Ketterle first 
{J*} ' showed had a two-step condensation: a GBEC into a band of states at a temperature T c and another 

"one-dimensional" transition at a lower temperature T\ into the ground state. In a thermodynamic 
limit in which the ratio of the dimensions of the anisotropic harmonic trap is kept fixed, T\ merges 
with the upper transition, which then becomes a normal BEC. However, in the thermodynamic 
i , limit of Beau and Zagrebnov, in which the ratio of the boundary lengths increases exponentially, T\ 

becomes fixed at the temperature of a true Type I phase transition. The effects of interactions on 
O |' GBEC are discussed and we show that there is evidence that Type III condensation may have been 

1 1 1 observed in the cigar trap. 
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^ \ I. INTRODUCTION 

Bose-Einstein condensation (BEC) has been intensely researched in recent years since the advent of laser and 
magnetic cooling of trapped alkaline atoms achieved that state [H-Q. In normal Bose-Einstein condensation (NBEC) 
a macroscopic number of particles populates a single quantum state (usually the ground state of the system) below a 
critical temperature. The state is often the lowest momentum state in a homogeneous system 

(IB (for 

an extensive 

review see Ref. 0), or it could be the lowest harmonic oscillator state in a trapped gas. Many authors have considered 

■ alternative possibilities: for example, a number of states may simultaneously be macroscopically occupied, or a band 
of states is macroscopically occupied while no single state has a macroscopic number of particles. Either of these 
last two cases is called a "generalized Bose-Einstein condensation" (GBEC). The terminology "generalized" was first 
used by Girardeau [§l-[l0| in considering a homogeneous interacting gas. (See also the work of Luban [ll|). An early 

— , description of the phenomenon was by Casimir |l2l) in treating a uniform gas in an elongated box (the "Casimir prism")- 
Rehr and Mermin [13| found that a rotating Bose gas had a GBEC. Other terminology has referred to "smeared" @ 
and "fragmented" condensations (l^ . [l5| . The most thorough analyses of GBEC have been done by Van den Berg and 
coworkers [l6l - [24j and more recently by Zagrebnov et al (25 33] . Ho and Yip [34[ claim that the spin-1 Bose gas is an 
example of fragmentation. 

However, there seem to be several schools of research on this topic who do not know of and do not quote the others' 
work. For example, Van Druten and Ketterle (3|| and more recently Beau and Zagrebnov [32] have theoretically 
discussed an example of GBEC. Ref. [32j has followed on the extensive work of the Van den Berg group and does 
not quote Ref. [35]. On the other hand, Ref. [35], and many references in the literature to this paper, do not refer to 
the previous papers of the Van den Berg or Zagrebnov schools. A frequently quoted paper by Nozieres (l5| presents 
a proof to show that fragmented BEC is ruled out for repulsive interacting systems; however, the proof holds only 
for one kind of GBEC; no proof or reference appears for the other kinds of GBEC. Thus we believe a pedagogic 
review of the subject is needed, to clarify the subject, make the details more widely known, and perhaps to stimulate 
experimental research to find cases of the phenomenon. 

Van den Berg et al [l9| have classified three general forms of BEC: There is condensation into 

I) a discrete number of quantum states each having macroscopic occupation (of order N, the number of particles) — 
e.g., NBEC has a single condensed state; 

II) a band of states with each state having macroscopic occupation. 

III) a band of states each state having only microscopic occupation (the occupation number rii of each such state 
has rii/N — >■ in the thermodynamic limit), but with the entire band having a macroscopic number of particles. 
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We will present examples of each of these situations in the present paper. All of our examples will be of ideal 
gases. The question then arises whether a GBEC can be maintained in the case of interacting systems. Nozieres (l5| 
showed that interactions would favor NBEC when there were repulsive interactions. Nevertheless there are several 
cases in the literature where interacting cases of GBEC have been given; we will discuss these later and how they 
avoid violating the Nozieres analysis. 

The treatments of Refs. [32] and [35] are interesting in that they bear on BEC in ultra-cold gases; they have 
theoretically discussed a case of two-step GBEC with a first transition into the lowest band of states in a "cigar"- 
shaped harmonic trap followed at a lower temperature by a condensation into the lowest single-particle state. The 
latter transition, into a single state, occurred at a transition temperature ~ 1/ ln(iV) , which would seem to disappear 
in the thermodynamic limit (TL). Each pair of authors considers an alternative TL in which the lower transition 
persists for large particle number. We will discuss these cases more fully below. Sonin (3(| was the first to note that 
there could be multiple BEC transitions in a parallelepiped geometry corresponding to a variety of way of taking the 
TL, although he did not discuss any associated GBEC. The original case of a GBEC associated with the double BEC 
transition was considered by Van den Berg and coworkers in a flat-plate geometry [H], [23| . There have been many 
experiments with gases in cigar traps; we comment on their relevance to the possible observation of GBEC. 

We begin by discussing the normal BEC (NBEC) for an isotropic harmonic potential followed by a sample of Type 
II condensation. Then we give a detailed discussion of the anisotropic harmonic "cigar" trap, for which transitions 
are possible at two different temperatures. As discussed above the nature of the lower transition depends crucially 
on how the TL is taken. Such a result shows a secondary purpose of our paper: to show there is more than one way 
to go to the large-scale limit; the density of states will depend on the relative values of the respective dimensions 
of an ansiotropic container. We treat the "Casimir prism" to illustrate GBEC with box boundary conditions in all 
directions. The effects of interactions on GBEC and the implications of recent experiments are discussed in Sec. I VII 
A more abstract analysis of the three kinds of GBEC is given in the Appendix. 



II. THE ISOTROPIC 3D OSCILLATOR 



To illustrate NBEC, we use the three-dimensional (3D) isotropic oscillator potential rather than the usual 3D 
homogeneous ideal gas, because it can show a very interesting GBEC when the potential is made anisotropic. The 
harmonic potential is taken to be 

V=^( X i+yi + z*), (1) 

where U is the potential strength. Note that the harmonic potential has been written so it has an evident length scale 
L and the frequency is 

m L 

where to is the particle mass. We insert this length since taking the thermodynamic limit in a harmonic potential 
involves increasing the number of particles while weakening the potential to maintain constant overall density p. 
Having this length scale is also necessary to see the relative sizes of the occupation numbers of various states [U EH . 
The harmonic length ao = \J h/(muj) is a measure of the rms deviation of a particle from the center of the well and 
is not the appropriate density length scale as we will see. The energy levels are 

£ P*PyPz = ?>U [Px+Py+Pz + T^J , (3) 

with pi taken over all positive integers. 

The total number of particles in the system is given in terms of the usual Bose distribution function as 

N= {cx P [/3(M^+P y +P 2 )-M)]-ir\ (4) 

P*,Py,Pz 

where /3 = l/ffc^T) with T the temperature, fc^ the Boltzmann constant, and /i the chemical potential. We will 
always consider the ground-state energy to be incorporated into the chemical potential ^ so that the ground state is 
effectively zero. We introduce a temperature parameter 

T _ h [ui . 

T ° = T~\ ' ( 5 ) 
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where a = 1/p 1 / 3 is an average interparticle separation. 

For high temperature T we can replace the sums by integrals: 
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du\ I du2 



e (u 1 +u 2 +u 3 +a) _ ^' 



where = T^api/TL and 

a = —ftp,. 

The integral is most easily done by expanding the integrand in inverse powers of the exponential to get 
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where the Bose integral [39] is 



F n (a) = £ 
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We see that the way to define a density parameter (or any thermodynamic variable) so that it is scale independent 
[37l |38| is in terms of a volume T 3 with the parameter given by 
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a 3 ~~ T 3 ' 



(10) 



As the temperature decreases, Eq. ([5]) can be satisfied by decreasing a until the quantity F^^a) has a maximum, 
^3(0) = C(3) = 1-202 where C( n ) is the Riemann zeta function. The temperature corresponding to this maximum is 
the critical temperature T c given by 



T c = 



To 



C(3)V3 



(11) 



Below the transition temperature the non-condensed particle number, N — No, with Nq the condensate number, is 
still given by the right side of Eq. (JBJ with a = : 



N-N, 
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from which we see that 
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Assume that all the condensed particles fall into the ground state (as usual); the occupation number corresponding 
to the ground state is 
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where O(x) means order of magnitude of x. Thus a is small, but not actually zero. The low excited states have 
occupation 



exp 
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YJ-VPx + Py + Pz) + a 



(15) 



Since a is so small we can drop that term in Eq. (P~5|) as long as some pi > 0. In that case the occupation number 
is 0(L/a) = 0(N 1 ^ 3 ), which is small relative to the ground-state occupation, O(N). Only the ground state is 
macroscopically occupied. All this is quite standard except here we have states of an oscillator potential rather than 
free particles in a box. The BEC considered here is Type I because only one state has macroscopic occupation. 

Normal BEC might also include cases where there are a discrete number of states each having macroscopic oc- 
cupation. For example, the well-know experiment (40| in which two condensates were brought together to form an 
interference pattern is such an example. One might also consider multiple condensates each having a different spin 
orientation [3J, [41| . Condensates trapped in multiple wells are sometimes said to be "fragmented" [42} . 



III. THE CHANNEL POTENTIAL 



In a Type II condensation we have a band of states each with macroscopic occupation. Here we discuss a rather 
peculiar case [l7j in which a spinless two-dimensional (2D) ideal gas has a harmonic potential in the z direction and 
is free in the x direction with periodic boundary conditions in that direction over a large distance L. The potential 
forms what might be called a channel or trough. The energy levels are 
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with the harmonic frequency given by Eq. ((5]); we use the same L parameter in both momentum and harmonic 
dimensions. The total number of particles in the system is given in terms of the usual Bose distribution function as 



N 
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exp 
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with s.. 



and 
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. We define two temperature parameters 
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(17) 
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where the density is p — 1/a 2 . These temperatures are simply natural units expressed in terms of the parameters of 
the problem. For simplicity we assume the two parameters are the same: T x = T z = T . For high temperature T we 
can change the sums to integrals: 



N = 



a 2 
N 
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T 
T 



dpz 
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(20) 



since the average density is p = N/L 2 . F 3 / 2 {a) has a maximum C(3/2) = 2.61 (as in the standard 3D homogeneous 
gas). The critical temperature is 



5 



For T < T c there is a condensate whose number Nq is given by 



N-N, 



n 



N 



N 



3/2 



V^F 3/2 (0) 
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or the condensate fraction fo = Nq/N is 



fo 



3/2' 



(22) 
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Assume that the ground state (that is, lowest k— state of the lowest harmonic level p z 
occupied with number 



0) is macroscopically 



1 
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(24) 



However, now the excited states cannot be neglected as in Sec. |TTJ Consider the entire band of k states corresponding 
to the lowest harmonic level p z = 0. The low excited states have occupation 
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where we have written 7 = aL 2 /a 2 . Obviously n^.o = 0(L 2 /a 2 ) for all k. On the other hand, we can easily see that 
the harmonic bands of states with p z > do not contribute macroscopically, being of order L/a <C N. We must sum 
up the contribution Nf, of the first harmonic band (p z — 0). We then have for the band number 



N =N b = 



OO 

E 

s x — — 00 
00 

E 



1 



exp [i(^ s | +7 )j 
1 



coth 




(26) 



The sum in the last line is exact; replacing the sum in the first line by an integral is not sufficiently accurate for all 
T (it sets the coth factor equal to one). Setting the temperature dependent result found in Eq. (f2"B")) to Nb/N from 
Eq. (|26p gives us an equation for 7, which can be used to evaluate the individual energy-state occupation fractions 
ris^fi/N. The results are shown in Fig. [TJ 

Near the condensation temperature the number of macroscopic states involved is clearly very large. Only near 
T = does the lowest state become dominant. All these occupation numbers are of order N. 



IV. THE CIGAR TRAP 



For another, and most detailed, example of GBEC we consider the anisotropic harmonic trap. A few years ago 
Van Druten and Ketterle (VDK) [35| suggested an unusual BEC that would take place in a cigar-shaped harmonic 
trap (one long dimension with a weak harmonic potential, plus two short dimensions with much stronger transverse 
potentials). More recently Beau and Zagrebnov (BZ) [32] independently treated the same problem. In each treatment 
there is a GBEC into a band of states, and a second transition at a lower temperature. In these treatments we 
learn that there are two ways to take the thermodynamic limit (TL); in one the transition disappears (that is, it 
is a pseudo-transition), while in the other it persists in the infinitely large system. This peculiarity illustrates some 
interesting unexpected features of the statistical mechanics of BEC. The approach used by BZ to fix this second 
transition was originally invented by Van den Berg and co-workers [2lL [23| for an anisotropic square well potential. 
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FIG. 1: The condensate distributions for the channel potential versus temperature. The solid line is the entire condensate 
distribution /o; the dashed line is the that of the lowest momentum state fo.o = no,o/N of the condensate band; the dotted 
lines are the distributions of the next higher momentum levels. 



The energy levels of the cigar trap are 



e PxP B P, = tWJ-(Px +Py) + tw\\Pz, 



(27) 



where u>i — y^U/m/Li corresponding to the two lengths L± and Lu, with Ln L±. The density parameter of the 
gas is then taken to be 



N 



LiL» 



(28) 



We again define a characteristic temperature To just as in Eq. (0. With this notation the particle number must 
satisfy 



"= E 



i 



Pm,Py,Pz eX P 

Change the sum to an integral for T large to give 
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Thus the chemical potential parameter a has to satisfy 

T 



F 3 {a) = 1. 



If T becomes too small, the equation can no longer be satisfied. The condensation temperature is 

T 



T r = 
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where the density p is given by Eq. 

The number of noncondensate particles now satisfies the relation 



V -Yn = V ( pj , (33) 



with No the condensate number or 



When the temperature is below T c the condensation is in a band of states, as we will show explicitly below. The band 
occupation number is 



^ = E — 1 — ■ ( 35 ) 

p 2= o exp 



T q p 2 
T L, 



If we change the sum to an integral we find 



TN f°° l 

N b = -— dw 36 

T ii Jo exp [w + aj - l 



in which we have made the substitution 



where 



L _ ( N \ l/3 _ iV 



,2/3 



^ = l^^T J ■ (3?) 

Note the assumption (to be corrected later) in this replacement of the sum by an integral, that the ground state 
occupation remains as small as that of the rest of the states. Carrying out the integration gives 

fb= — = -—\n(l-e- a ). (38) 

If indeed the ground-state occupation is negligible, then the band occupation is the whole condensate, fb = /o and, 
since a is very small below the transition temperature, we have ln(l — e~ Q ) In a; we can solve for a to find 



-foQK 



(39) 



The band states are analogous to a one-dimensional (ID) gas of Bose particles in the long harmonic trap. Does 
this ID gas also have its own condensation into its lowest-energy state, which is, of course, the overall ground state? 
Subsequent treatment now depends on the boundary conditions used in taking the TL. 



A. Standard boundary conditions 

Under normal thermodynamic limit conditions, the ratio 

A = -I (40) 

is held constant, while N, Lii and are increased. In the present TL, K = (N/ A) 2//3 (Eq. (f3"T)0 would be increasingly 

large as N increases (~ N 2 ^) and, for T < T c , a would be increasingly small going asa~ e -JV " /3 from Eq. (|3^|) . 
The occupation numbers of the low-lying band levels are 

rp jy T U 
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The first term in the denominator being of order iV -1 / 3 dominates a giving uqo Pz — 0(N 1 ^ 3 ) and 71qq Pz /N — > in 
the TL (microscopic occupation); the condensation is then of Type III as listed in Sec. |TJ The p z = term is the 
exception; when iiqqo/N is microscopic, the condensation remains of Type III. However, the ground state occupation 
can become O(N) when 



n 000 



- = e^ K = cN, 



(42) 



where c < 1. This occurs at a temperature T\ given by 

Ti 



To, 



K 



/o(Ti) ~ In cN ' 

or, assuming N large and T± is sufficiently smaller than T c that fo ~ 1, we have 



ft 



E/ N 



m LuhiN 



(43) 



(44) 



which might be considered the condensation temperature of the ID band. In the TL, T\ would go to zero as 1/lniV 
if the ID density N/L\\ were held constant. Such a behavior is indeed characteristic of a ID pseudo-phase-transition 
[38l . |43| , which does not exist except in a finite system. While it is possible to consider a 3D TL such that L± is held 
constant while N and L« increase, it is p = N/(L^_L^\) and A that are held constant in the standard 3D TL while the 
two lengths increase. In Eq. (|4"3")l . when the right side becomes very large, ~ N 2 ^ 3 /\nN, the only way the equation 
can be satisfied is to have T\ ~ T c so that the denominator on the left, fo(T\), almost vanishes (See Eq. ([34])). Thus 
in the TL, T\ — > T c and the "ID" transition coincides with the condensation into the band of states, that is, the band 
then collapses to just the ground state and the system reverts to the normal BEC rather than a GBEC. Such a result 
is rather surprising and has not been noted before in the literature. See below for a numerical example. 

When the ground state occupation becomes large, the analysis of Eq. (|36p is no longer valid, and we need to be 
a bit more careful to find an expression for the ground-state occupation number. In the sum of Eq. (|35l) we should 
split off the ground state term n„ to write 



iV = n q + N b 



p^o exp 



loKEk +a > 

T N 



- 1 



(45) 



where we have redefined the summation index p' z = p z — 1 and let 
The ground-state occupation fraction is 



(46) 



/s iV(e a -l)' 



Replace the sum by an integral as before to give 

A 1 1 
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In 
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(47) 



(48) 



Since the total condensate fo is given by Eq. (|34j) . this is a transcendental equation for a. 

A formula to solve self-consistently for f g for all situations is gotten from Eq. (|4"5)) by solving Eq. (|4"T|) for e a in 
terms of f g : 



fg = f + 



T 



In 



1 — e nt 



1 



(49) 



In the TL the factor in front of the logarithm in Eq. (l49l) dominates the behavior and the second term goes to zero 
so f g —> fo as we stated above. 
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B. Numerics 

Generally we have 

VDK considered a set of parameters with N = 10 6 and A = 5.6 x 10 4 . In the this case we have 

T X 2/3 , 1q6 \ 2/3 



The numerical result in Eq. (f52"j) was gotten by iteration; if we set fa — 1 in the T\ formula we get T\/T c = 0.52; 
putting that back on the right in /o(Ti) gives a new value of Ti, etc. In Fig. [5] we plot the results for fg and f g , 
where we have solved Eq. (|4U1) by simple iteration. 
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FIG. 2: The condensate distributions in the standard boundary limiting case (the ratio of boundary lengths A = Lu/LxiB a 
fixed constant) and for N = 10 6 , A = 5.6 x 10 4 . The solid curve represents the full condensate fo (the fraction of particles in 
the lowest band of states), while the dotted curve is the ground-state pseudo-condensate f g (the fraction in the lowest state of 
the lowest band). 



In the case where we keep the same length ratio, but increase the number of particles to N = 10 8 we have 



K = 
T, 



10 



8 \ 2/3 



5.6 x 10 4 

/ (Tj)irC(3) 1/3 
In TV 



= 147, 
= 0.961 



(53) 
(54) 



Fig. [3] shows the plots of the two condensate fractions for N — 10 s . The ground-state transition almost coincides 
with the band condensate as expected from the analytic argument. 

Ti approaches T c because in Eq. (1431) the right side becomes very large with K ~ TV 2 / 3 — > oo, requiring increasingly 
large l//o(Ti) on the left as explained above. The question then arises whether there is a TL such that the right side 
does not become large, that is, in which K ~ In N . In the next section we will find such a case. 



C. Exponential boundary conditions 



Here we consider how the above discussion is changed if we now take the length ratio to obey 



(55) 
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FIG. 3: The condensate distributions in the standard boundary limiting case (the ratio of boundary lengths A = L^/L ±is a 
fixed constant) and for N = 10 s , A = 5.6 x 10 4 . The solid curve represents the full condensate /o, while the dotted curve is 
the ground-state pseudo-condensate f g . 



a boundary condition that certainly forces the system to become more ID as L± increases. Similar exponential 
boundary relations were first proposed by Van den Berg and co-workers (HHHI for an anisotropic square well potential. 
But these explicit conditions for the anisotro pic harmonic oscillator, in which the length factor in the exponent is 
squared were given by Beau and Zagrebnov |32l . [33j and we designate them as BZ conditions. That they fit the 
requirement that K ~ In N will be verified below. It becomes useful to use a unitless notation here. Let 



t± = P 1/3 L ± , 
9 

,2/3 



7 = 4m- ( 56 ) 



Then, from the density relation Eq. ([28]) . we have 



K = r ' 2 



7 - ^ (.7) 

The VDK parameters TV = 10 6 and A = 5.6 x 10 4 can be considered a special case (one particular N value) of a BZ 
set. Then N and A give the corresponding 7 value. From the above formulas we have £± — 2.61 and 7 = 1.60, which 
is the 7 value we will use later for larger N values at the same density. For the original VDK parameters, the curves 
for /o and f g are obviously identical to those shown in Fig. (5J however, with the condition Eq. (|55p and a larger TV 
value, f g will no longer merge with the fo as we will see. We consider a much larger N value and solve Eq. P^)) for 
f g . We show the comparison of the f g curves for N = 10 6 and N = 10 16 in Fig. 0] We can understand the behavior 
of the large N curve by approximating Eq. (|49[) : For large N (any we have used here) we can expand the exponential 
and the fraction inside the logarithm to find 



T 

fg = J0+ In 
1 K 



TqK , 1 



NT Nfg 



/o--^ln[iV], (58) 

IqK 



where the last approximation holds because N is much larger than the other factors in the logarithm. The parameter 
K = i\ can be expressed approximately as well. From Eq. (|57[) we write 



j£\ = InN - -ln£ 2 x . (59) 



By iterating this formula once, we see that the second term is much smaller than the first and so 



(60) 
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FIG. 4: The ground-state condensate distributions f g in the BZ boundary limiting case (the ratio of boundary lengths L\\/Lx 
e 9ij - is not a fixed constant in the thermodynamic limit). The dotted curve is for N = 10 6 while the solid curve is for N = 10 1 



Putting this in Eq. 



gives 



f g (T) = f (T) 



T 7 



T C C(3)V3 



(61) 



We plot this result in comparison with the N = 10 16 result in Fig. 
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FIG. 5: The ground-state condensate distributions f g in the BZ boundary limiting case for N = 10 16 (solid curve) and the 
thermodynamic limiting value, Eq. (|6ip (dotted curve). 



The lower transition temperature is gotten by setting the left side of Eq. (|61|) to zero: 



Ti 



T c /o(T 1 )C(3) 1 / 3 
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(62) 



which must be solved self-consistently. We find Ti = 0.552T C . Clearly this approximation is very good and the system 
has two distinct phase transitions as BZ have claimed. Perhaps we should not be too surprised that this was possible. 
In the VDK case having T\ « 0.5T C was arranged by a judicious choice of A for N = 10 6 . The question then is 
whether it is possible, for a larger value of N, to find a A such that Ti and p are both unchanged. Such a value is 
given by Eq. (|55l) . VDK also numerically treated a case in which T\ was fixed while N increased, although they did 
not specify how this was done; they apparently used a form equivalent to Eq. (|55[) . 

There is a case of box boundary conditions (square- well potential) that is isomorphic to the harmonic system treated 
in this section. This is the flat plate geometry in which two large square plates of length L on a side are separated by 
a much smaller distance D. This problem was treated by Van den Berg and coworkers [2l|,|23j in the 1980's. A third 
analogous case involves having two free dimensions and one harmonic dimension [32] . 

Sonin [36] gave an analysis of multi-step quasi-condensations in finite systems for anisotropic free-particle boundary 
conditions and showed that a second transition could be preserved in the thermodynamic limit by an appropriate 
thermal limiting procedure. Deng (47| and Shiokawa (48| treated multi-step transitions in finite anisotropic harmonic 
potentials; the extra transitions disappear in the normal TL. 
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V. GBEC IN A BOX: THE CASIMIR PRISM 

We consider one more Type III GBEC where there is again a macroscopic condensation into a band number 
of states while the occupation of any single quantum state remains microscopic. With th is g eometry there is no 
second transition at a lower temperature. The case considered here was treated by Casimir [12j and later by others 
[TgL l23l l30l . HH]. It was one of the first known theoretical cases of GBEC. The geometry is shown in Fig. [5] The 
length L is much larger than the side D of the square cross section. 




FIG. 6: Casimir Prism. A box with one side L much greater than the square end face dimension D. 
The free-particle states in this system are given by 



2m 



(4 + Sy) _ 

D 2 + L 2 



where the Si are again positive or negative integers. The total number of particles is a sum 

1 



Above any transition we can change the sum to an integral: 



L 



I)-m]} 



- 1 



N 



00 poo 

ds x I dsy 

— 00 J — 00 



ds 7 



LD 2 {2m\ 3/2 f°° , f°° , 

— — / dui / du 2 



1 



/1 3 V /? 

^3/2(0;) has a maximum of C(3/2) at a = so the transition temperature is given by 



2/3 



27TTOfc B VC(3/2) ; 

For T < T c , Eq. ([55)1 is no longer valid, and below T c the condensed particle number Nq satisfies 



N Q = N 



3/2 



(63) 



(64) 



(65) 



(66) 



(67) 



13 



Because of the anisotropy of the boundary conditions, a band of states with s x = 0, s y — fills at the transition. To 
see how this occurs, we examine the density of low excited states. There are various ways to take the thermodynamic 
limit, but here we hold D constant while we take L —> oo. (The arguments presented here also hold if we let L, D — >• oo 
with L approaching infinity faster than D 2 . Compare with Appendix A. 3.) Let us first hypothesize that the ground 
state is macroscopically occupied so that a ~ 0{1/{D 2 L)) (which we will see is incorrect); the occupation of the low 
states would then be 



1 



h 2 /3 
Ira 



' («;+«;) 

D 2 



L- 



+ a 



(68) 



Under our hypothesis, the first term would dominate the other two, which could be dropped and the state density 
would be 0(1) unless s x = s y = 0. In that case, the s z term would still be negligible and we would find n fl,s z ~ O(N) 
for all values of s z . Clearly we need to sum up the entire band of s z states to find the true behavior. The number of 
particles in the band is 



s z exp 



/3h 2 s 2 
2m T 1 



2m 
2nm 



1/2 



L <-oo dU 'exp{w 2 +a}-l 



1/2 



LF 1/2 (a). 



(69) 



Below the transition temperature, a is small and Nb equals the condensate number of particles A^o- Further, since 
f39l | for small a, F 1 / 2 {a) w \/n/a, one gets 



or 



2n 2 m 
(3h 2 



1/2 



L 

Ja 



(70) 



2n 2 m ( L 



ph 2 \N 
2n 2 m 1 



(71) 



where p — Nq / LD 2 . This result is quite different from what it would have been if the ground state were macroscop- 
ically occupied as in our initial hypothesis of a ~ 1/(D 2 L). 

The condensate number is macroscopic, but the numbers in the single-particle states are not. The density corre- 
sponding to low quantum numbers {0, 0, s z } is 



1 



Po.o,. 



LD 2 



h?P sj 7 (T) 
2m TP ' ~LW 



h 2 /3 2D 2 , ~«L 
2m Z L "r" D 2 



(72) 



where 7 = aD A . The term in 7 always dominates so the occupation of each of these states is microscopic and 
approaches zero as 1/L, including the ground state with s z = 0. Of course, at T — the ground state must finally 
have all the particles in it, but the temperature at which this occurs can be estimated by using a from Eq. ()71j) in 
Po.0,0 ~ l/(Na) to see that the onset temperature for macroscopic occupation of the ground state is T ~ T' c p 1 / 3 D 2 /L, 
which is zero in the TL or very small in a real experiment. 

There is an alternative way to take the limit in which all three dimensions of the box become infinite. This way 
can also be used to distinguish the three kinds of GBEC. We discuss this in the Appendix. 
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VI. THE EFFECTS OF INTERACTIONS ON GBEC 
A. GBEC and interactions in the literature 

The purpose of this paper has been to outline the possible kinds of GBEC by using ideal gases. It is not our purpose 
to make a complete analysis of the existence of GBEC with arbitrary interactions. Nevertheless it makes sense to 
ask whether GBEC would disappear with interactions. Nozieres [15| has shown, in the Hartree-Fock approximation 
in a homogeneoous scalar Bose gas that repulsive exchange interactions favor condensation into a single state. He 
ignored the case of attractive interactions, because they would cause the homogeneous system to collapse. More 
importantly, he also assumed that each condensate state is macroscopically occupied, which applies only to Type I or 
II condensation. However, in a trap a small degree of attraction does not necessarily lead to collapse [44|; the kinetic 
energy stablizes the system. The literature also contains examples of GBEC's in interacting systems [25I4291 . l45l |4^| . 
involving repulsions and Type III band occupation. Girardeau [10] considered an attractive interaction in a uniform 
system, which showed macroscopic occupation of each condensate state, but he did not take into account the possible 
collapse of the system. 

A common type of interacting model showing GBEC has diagonal interactions. These interactions are a function 
only of the number of particles Nje in the fcth momentum state. Then the Hamiltonian is a function of a set of 
mutually commuting operators with a particularly simple spectrum. (See for example, Ref. [25j and references 
therein.) While one finds Type III GBEC in the interacting case, interactions also introduce yet another type of 
BEC called non-convential or dynamic condensation. The conventianal condensation occurs when there is a kind of 
saturation: the total particle number becomes larger than some critical value as in the NBEC or even Type III. A 
dynamic condensation occurs only when induced by attractive interactions. 

B. The Hohenberg theorem 

We treated only one case of a 2D system, that of Sec. IIIII where we found a Type II condensation, that is, 
having a band of macroscopically occupied states. However, the Hohenberg theorem [49] states that no macroscopic 
condensation can occur into a zero-momentum state in two dimensions in the thermodynamic limit. Since the theorem 
refers to condensation into a momentum state it might not seem to apply to the condensation discussed in Sec. IIIII 
since one dimension, at least, involves a harmonic potential. However, a theorem developed by Chester [50] based on 
work by Penrose and Onsager [5l| states that there can be no condensation into any state unless there is one into 
the k = state. The loophole relative to the trapped gas is the condition assumed by the Chester derivation that 
the density be finite everywhere in the thermodynamic limit. As we show below, the gas studied in Sec. IIIII has a 
divergent density at the origin in the TL. If, however, the system has repulsive interactions, no such divergence would 
be allowed and then the Hohenberg theorem would apply and the transition studied in that section would disappear. 
A similar situation has been discussed for a rotating Bose gas [l3[ and with the 2D completely trapped gas (38j. 

To see the divergence at the origin for the system in Sec. IIIII we consider just the ground state contribution to the 
density. If ipo{z) is the ground-state harmonic wave function, then 

fn\ 1 I / fn\|2 /mw\i/2 n 00 

Po(0) = — |Vo(0)| n 00 =l^—) — 

since noo = O(N) and L = \/N/p. In the TL, the density diverges, and the Chester theorem does not apply. 

The next question then is whether a Hohenberg-like theorem forbids a 2D transition of Type III, in which a band 
of microscopically occupied states condenses. This subject has been addressed and the usual derivation was 
shown not to forbid such a transition. However, the derivation in this reference does not tell us whether some other 
theoretical approach might not reveal such an alternative theorem forbidding the 2D transition. The question of 
whether a Type III GBEC is possible in 2D seems still an open question. Recently an analysis [53[ of an interacting 
2D trapped gas showed no fragmentation. In ID an analysis similar to that in Ref. [52| shows that no GBEC can 
occur in the standard type of thermodynamic limit. 
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C. The 3D transition to ID 

One of our prime examples of GBEC in Sec. IIVI involved the crossover between a 3D and a ID gas in a cigar trap. 
The literature on bosons in cigar traps is much too large to summarize here and no definitive GBEC has yet been 
seen experimentally. A good review is the paper of Bouchoule et al [54[ and a recent relevant experiment is described 
by Armijo et al (55|. Actual experiments here are on finite systems, of course, and sharp transitions are not observed. 
One theoretical advantage of ID is that there is an exact analytic solution due to Lieb and Liniger (5(| for bosons 
interacting by a S- function interatom potential, with thermodynamics by Yang and Yang (57| . This theory has been 
used extensively in analyses of the experiments. Forrester et al [58| have been able to show that, for the ID Bose 
gas with infinite (5-function interaction, whether homogeneous or harmonically trapped, all the eigenvalues of the 
one-body density matrix at T = are of order y/N. Such a result would certainly correspond to a GBEC if we could 
state that the system was a quantum fluid. But this impenetrable gas would seem more like a solid than a fluid! 

Armijo et al [55| map out the dimensional crossover from a 3D gas to a ID gas in a cigar trap. They observe 
a transition, which scales experimentally exactly according to the description in which "atoms accumulate in the 
transverse ground state, although no single quantum state is macroscopically occupied." This is precisely what we 
describe as GBEC in Sec. IIVI Theoretically no lower sharp second transition to a true BEC, in which only the ground 
state is occupied is expected, but rather a "quasicondensate" is formed. The formation of the quasicondensate is driven 
by interactions, which surpress density fluctuations while the phase still fluctuates [59j . The observations conform to 
this description. 

In this regard the path-integral-Monte-Carlo (PIMC) calculation of Nho and Blume [6(| on the 1D-3D crossover 
is particularly relevant. They compute the superfluid component in the gas. With the PIMC approach it is difficult 
to compute the occupation numbers for the interacting gas. However, using those of the noninteracting gas, they 
find that the superfluid component tracks the occupation of the entire band corresponding to the lowest transverse 
harmonic state much more accurately than the occupation of the lowest level in that band. 

A further indication is the work of Witkowska et al [6l| in which the authors use a "classical field approximation" to 
study the evaporative cooling dynamics of a trapped interacting ID Bose system. They compute the eigenvalues of the 
one-body density matrix and find, in an intermediate temperature range, that the lowest four states have occupations 
of ~ 10% or more. As the temperature goes even lower only the ground state remains occupied. This looks much like 
a GBEC or at least a "quasi-generalized-condensate." 

D. Spin-1 Bose gas 

Ho and Yip [34j have described a spin-1 Bose gas with antiferromagnetic spin-spin interaction, which has a frag- 
mented condensate ground state. Here angular momentum conservation prevents spin flips between +1 and —1 states. 
The ground state is analogous to having a scaler Bose gas in a double- well potential. A magnetic field gradient would 
allow such spin transitions and results in a non-fragmented condensate state. 

VII. CONCLUSION 

Our purpose in this paper has been to give a tutorial on generalized Bose-Einstein condensation to clarify what 
seems to be confusing literature on the subject. The research in the field has gone in directions with some researchers 
having been quite unaware of the results found by others. Proofs that GBEC cannot exist, which, to some readers 
seemed general, simply do not apply to other forms of the phenomenon of which the authors were not even aware. 

We have seen that there are three types of GBEC and have given examples of each. We have examined a form 
particular relevant to recent experiments, that in the cigar trap, developed independently in Refs. [32| and [35j and 
have seen how the properties in the thermodynamic limit can be completely different depending on how that limit 
is taken; under a rather pecular limit there can be a two-stage condensation in the ideal gas. It is possible that 
experiments in cigar traps have already observed the upper transition, although it is unlikely that there could be any 
sharp lower transition because of interactions and because no thermodynamic limit is actually taken. 

Our hope is that this paper might stimulate further research in this area, especially in interacting systems, so that 
ultimate experimental verification of the existence of the phenomen might be observed. 
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Appendix: Alternative set of boundary conditions 

Van den Berg and co-workers [23j have used a uniform approach for general boundary conditions that allows one to 
distinguish three geometries under which the three kinds of BEC occur. However, these geometries are not in every 
case identical to the situations we have discussed above as we will see. Consider a general box of sides Li,L 2 , and 
L3; all three of these dimensions will be taken to infinity in the thermodynamic limit. The TL involves defining an 
arbitrary length parameter of, say, atomic size (we choose the interparticle separation a) and a unitless parameter H 
that will be taken to infinity to establish the thermodynamic limit. We define 



Li = aH Vi , 

where the Vi are fixed parameters such that the volume V of the box is linear in H : 

V = LiL 2 L 3 = a 3 H, 

so that 

V\ + l>2 + V?, = 1. 

The parameters are arranged according to 

v\ > v 2 > V3 > 0, 



(74) 
(75) 
(76) 
(77) 



and the three types of GBEC are categorized according to whether v\ is smaller than 1/2, equal to 1/2, or greater 
than 1/2. 

We will assume we are below the 3D transition temperature. The density of the lowest states is given by 



1 



V 



a 3 H |exp 



Ph 2 
2ma 2 



h 2 /3 



2 m 



where 



V 

7 = a-3. 



In every case the condensate density satisfies the usual 3D behavior, 



Pa= P 



3/2 



(78) 



(79) 



(80) 



as shown in, say, Sec. [V] 



A.l Type I. u ± < 1/2 



The NBEC in a cubic box corresponds to V\ = U2 = v§ = 1/3. More generally we see that for v\ < 1/2 every 
1 — 2ui > 0. Assume the smallest of these is 1 — 2v\. We hypothesize that the ground state number is macroscopic 
so its density satisfies a 3 p oo = 1/7 = 0(1). Then in Eq. (|78|l the terms in s dominate over the 7 term and an 
excited-state density is - l/i/ 1 " 2 " 1 (if s 1 = 0, then p 0jS2 , S3 is 0(l/i7 1 ~ 2l/2 ) or 0{l/ H 1 ^ 3 ), which is even smaller). 
In the thermodynamic limit we have H 00 and the excited state densities vanish. Thus only the ground state is 
occupied and our hypothesis is verified. 



A. 2 Type II. m = l/2 



In this case 1 — 2vi = 0. We again suppose that 7 = 0(1), so that the dominant term in the square bracket of 
Eq. ([75)1 is either that in s 2 or S3 of order H 1-2 ^ 2 3 so that the corresponding density vanishes as H — > 00. However, 
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if S2 = S3 = 0, then the term in si is the same order as 7 , so that there is macroscopic occupation of each state, and 
we need to consider the whole S2 = S3 = band of states. This band has density 



1 



Pband 



h 2 2 1 



(81) 



In the present geometry it is no longer valid to replace the sum in Eq. (|8ip by an integral since both terms in the 
denominator are of the same order of magnitude. Here we must solve Eq. (|8ip for 7. The sum can be done analytically 
in terms of a hyperbolic cotangent (Sec. IIII|) giving a transcendental equation for 7. 



A. 3 Type III. u x > 1/2 



Now we have 1 — 1v\ < while the other two such factors are positive. The state density can be rewritten 

1 

Psi,S2,S3 ~ 



(82) 



where each of the exponents is positive. In Type III all states are microscopically occupied so we should try assuming 
that 7 = H v where r] is a positive number less than 1 — 2v?,. In that case the order of magnitude of the densities with 
si is either l/H 1 ^ 2 " 3 or l/iJ 1-2 " 2 , the latter case occurring only if S3 = 0. These densities vanish faster than 
the densities of the band of states with S2 = S3 = 0, si ^ 0; the latter are all of the same order of magnitude l/a 3 7 
as the ground state and we need to sum these band states to get the entire condensate: 



Po 



E 



1 



h 2 /3 2 , 



1 (2ma 2 H 2v ^ 



h 2 p 



1/2 



7T 

coth 



/2ma 2 i? 2 ^- 1 7 



V7 



h 2 /3 



which yields 



7 = 0(H 



2i>i-l\ 



(83) 



(84) 



so 77 = 2vi — 1. In the limit of H — > 00, p Sl .o,o = 0(i7 1 2l/1 ) —> so each condensate state is microscopically occupied 
as assumed. 

' h 2 0k 2 



In Eq. 



the distribution is 



(^+7) 



, which will become negligibly small for momenta at k > fco with 

h 2 pk 2 



2m 



(85) 



where C is, say, 10 4 . The cutoff momentum is 



fc ~ 1/H 



l-Vl 







in the TL; the condensate bandwidth in momentum space is vanishingly small. Indeed Girardeau defines GBEC in 
the following way for a homogeneous system He writes 



f = lim lim — nk, 



(86) 



fe<fco 



defining "/ as the fraction of the total number of particles with momenta less than any macroscopic momentum" and 
where nu is the number of particles in momentum state k. 

However, the quantum number corresponding to this cutoff momentum is 



s n 



afco 
"2T 



H V1 -H 2 " 1 - 1 -> 



(87) 
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since v\ > 1/2. So while the condensate band has all its significant momenta microscopic in the TL, the number of 
such states in the condensate is infinite! 
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